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Похибка трансформацій аномалій МГПСТ, як гармонічних функцій, залежить від кривизни еквіпотенціаль- 
них поверхонь поля, амплітуди аномалій, міри області трансформації, 1 сягає, залежно від виду трансформації, 
недопустимих величин. Інтерпретація гравіаномалій як гармонічних функцій неефективна для регіональних 
структур. 
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Гравіметрія, як і ряд інших наук про Землю, підійшли в своєму розвиткові до етапу, на якому 
отримання нових узагальнених знань про глибинну будову земної кори і мантії неможливе, як 
відзначають автори |Страхов и др., 2001; Черньтгй, 1992; Интерпретация..., 1992) без постановки 
глобальних експериментів, що вимагає обгрунтування застосовності наявних методів інтерпретації 
потенціальних полів до глобальних масштабів. Формалізація моделей середовища та апріорної ін- 
формації відбувається на основі апроксимаційного підходу |Страхов и др., 2001). Особливої значу- 
щості набувають проблеми повного видобутку корисної інформації за рахунок редукцій, транс- 
формацій 1 комплексування методів |Гравиразведка..., 1990) та побудови адекватної глобальної 
нормальної моделі Землі, щодо якої здійснюються ці трансформації. Постає питання роздільної 
здатності методу. 

Актуальність повного видобутку корисної інформації з даних гравіметрії засвідчує чимало пу- 
блікацій |Оганесян, 2004; Балк, 1989; Страхов, 1990; Старостенко 1973; Дубовенко, 2009а| (як 1 на- 
ведена в них бібліографія). В руслі вивчення цієї проблеми формується нове відгалуження в теорії 
інтерпретації даних гравіметрії - теорія відновлення модуля градієнта потенціалу сили тяжіння 
(МГПСТ) (|Черньгй, 1992|; в її рамках досліджена нова некласична задача відновлення МГПСТ 
| Дубовенко, 2009а; Дубовенко, 20096; Дубовенко, 2010). 

19. Актуальність проблеми. В процесі інтерпретації гравіметричних спостережень чільне 
місце посідає редукція сили тяжіння за допомогою методів аналітичного продовження, серед яких 
особливо зауважимо задачу відновлення потенціалу в заданій області. Темі редукції сили тяжіння 
присвячена численна бібліографія |Гравиразведка..., 1990; Старостенко, 1978; Дубовенко, 2009а; 
Дубовенко, 20096; Дубовенко, 2010; Черньгй, 1969; Алексидзе, 1965; Черньй, Якимчик, 1999; 
Страхов, 1964|, але досі вона вирішується наближено через неточність постановки класичних за- 
дач теорії потенціалу. Дійсно, при розв'язанні зовнішньої задачі Неймана для рівняння Лапласа, як 
і редукції за методом Стокса-Молоденського, маємо точний диференціальний оператор для потен- 
ціалу сили тяжіння, але наближені граничні умови: земна поверхня не є еквіпотенціальною, а сила 
тяжіння не є лінійною комбінацією гармонічної функції та її нормальної похідної. В редукції за 
допомогою розв'язання задачі Діріхле для рівняння Лапласа навпаки, маємо наближене диферен- 
ціальне рівняння, оскільки величина РЕЗ; негармонічна | Алексидзе, 1965|, зате точні граничні 


умови. Першому методу надають перевагу в геодезії, а другий є традиційним для гравіметрії. Така 
підміна виявляє недоліки класичної редукції, якими можна нехтувати в малих масштабах, але які 
стають відчутними уже в регіональних побудовах. Нагадаємо основні вади редукцій: 

1. розрахунок нормального поля Землі 1 поправок за висоту точки спостереження виходять з різ- 
них моделей нормальної Землі; 

2. слабко враховано розподіл густини в поправці за проміжний шар 1 топопоправці (за рельєф); 

3. аналітичне продовження аномалій сили тяжіння як гармонічних функцій дійсне лише для об- 
ластей малої міри |Черньгй, Якимчик, 1999|; 

4. неврахування різної направленості векторів нормальної УС) і реальної «(х) сили тяжіння 
провокує, як покажемо нижче, похибку майже в 200 мГал при максимальному розходженні векто- 


рів до 407, а всередині такої похибки може сховатись реальна аномалія від глибинних джерел. 

Ці фактори обмежують сферу застосування наближених по-суті класичних методів редукції 
областями малої міри. Іншими словами, традиційні методи гравіметрії, основані на моделях плос- 
кої Землі і гармонічних аномалій сили тяжіння, заданих на необмежених множинах, не мають на- 
лежної роздільної здатності для інтерпретації гравіаномалій від глибинних аномалієтвірних дже- 
рел, незважаючи на їхнє поширення в геофізичній практиці. Мода на подібні глибинні густинні 
побудови, як зазначає |Черньгтй, Якимчик, 1999|, веде до неоправданої дискредитації методів аналі- 
тичного продовження в очах кінцевих користувачів через їх уживання в неналежних геологічних 
ситуаціях (зачаровані магією сучасної наукової графіки, інтерпретатори інколи забувають, що в 
основу тлумачення закладені прості моделі середовища з усіма властивими їм припущеннями й 
обмеженнями). Внаслідок значущості похибок перерахунку аномалій класичними методами для 
вивчення глибин Землі, очевидно, необхідно задіяти повнішу інформацію про поведінку аномалій 
сили тяжіння. 

Вихід з цього тупика неможливий через просте вдосконалення алгоритмів продовження, а 
пролягає через указання меж застосовності стандартних методів |Серкеров, 1990|, чисельну оцінку 
похибок. Це дозволить розділити ,,сфери впливу" класичної теорії потенціалу 1 ,/глобальних" мето- 
дів. 

29. Витоки проблеми. Класичну методику відновлення потенціалу закладено в | Юнг, 1936|, 
де запропоновано вважати прирости МГІПСТ в локальній області , вертикальною складовою сили 
тяжіння, і для наближеного аналітичного продовження сили тяжіння в малій області переформу- 
лював задачу з правою частиною у вигляді наближених значень сили тяжіння в задачу Неймана 
для рівняння Лапласа, що дозволило задіяти для аналітичного продовження апарат теорії гармо- 
нічних функцій. З кінця 1950-х рр. за допомогою аналітичного продовження досі вирішують задачі 
локалізації аномальних джерел далеко не завжди в локальних областях, про пагубність чого вже 
згадувалось. 

Перелом міг настати з виходом праці | Алексидзе, 1965, в якій вперше зауважено, що аномалії 
сили тяжіння загалом не є гармонічними функціями, якими їх в першому наближенні вважав Юнг, 
але праця пройшла повз увагу геофізиків. Основні труднощі редукування полягають в тому, що 
для повного опису значень поля сили тяжіння неможливо отримати лінійне диференціальне рів- 
няння, якщо невідома форма й густина збурювальних тіл, тому доводиться відшукувати деякі на- 
ближення збурень гравіполя. Редукцію саме в такій постановці досліджував Алексідзе, відшукую- 
чи гармонічне наближення реального розподілу сили тяжіння в деякій обмеженій замкнутій обла- 
сті, розміри якої узгоджені з точністю вхідних даних. У побудові цієї локальної області вжито 
властивість гармонічності похибки апроксимації значень сили тяжіння. Звісно, точність апрокси- 
мації обернено пропорційна розмірам області, в якій апроксимується сила тяжіння. Подробиці 
опису гравіаномалій в глобальній області дає |Чернь(й, 1982). 

У силу сказаного доцільно переглянути методику відновлення потенціалу в глобальній облас- 
ті за його слідами (чи його похідних) на межі області (задачі Діріхле, Неймана чи Коші |Левенков, 
19871). Дійсно, гармонічний потенціал однозначно відновлюється в глобальній області при точних 
граничних умовах, як розв'язок відповідної граничної задачі для рівняння Лапласа |Сретенский, 
1946), але саме точні умови недоступні в силу методичних чи економічних міркувань". Лишається 
лазівка у вигляді розв'язання задачі Пуанкаре (про косу похідну), але для її вирішення не вистачає 
граничних даних про орієнтацію вектора сили тяжіння щодо земної поверхні. 

39. Розбіжність векторів. Напруженість поля сили тяжіння має вигляд 
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| Частинні похідні потенціалу за будь-яким напрямом співпадають з проекціями аномалій за цим напрямом, тобто доли ОМ, дез 
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спрямоване по нормалі у до рівневої поверхні еліпсоїда - генератора аномального поля у. На площах діаметром в кількадесят км 

зміною напряму у можна знехтувати і тоді аномалія бе. буде наближено гармонічною функцією, оскільки, через неперервність 


других похідних потенціалу У/, маємо Або ЗА У д ЛУ 2 0» ЩО ДОЗВОЛЯЄ відновлювати потенціал, як гармонічну функцію. 
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Задіяти виміряні значення сили тяжіння як точні граничні умови неможливо через їх негармонічність, негладкість та обтяження 
різноманітними похибками. 
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де чн соб(п, х; ), і - 1,3 - напрямні косинуси одиничного вектора п (х) внутрішньої нормалі до 
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Рис. 1. Приріст кутів між еквіпотенціальними поверхнями та фізичною поверхнею Землі в різних точках профілю. 


еквіпотенціальної поверхні «АМУ(х): У(у)- Сх у точці х. Очевидна ортогональність напруже- 
ності поля до поверхні (2), тобто вимірюється насправді приріст градієнта потенціалу, а не 
один з його складників, як це зазвичай уявляють |Веселов, Сагитов, 1965). 

Пояснімо це твердження детальніше. Під час вимірювань гравіметр встановлюють в кожній 
точці х, профіля бу за рівнем на еквіпотенціальній поверхні (2) - Сх, яка перетинає бу у цій 
точці (рис. 1). Просторову орієнтацію гравіметра при вимірюванні не фіксують, тому нахили 
приладу при зміщенні по профілю ду лишаються невідомими. Різниця цих нахилів зумовлена 
різною кривизною кожної з непаралельних одна одній еквіпотенціальних поверхонь (2) - Сх, 
Сх, зе Сх, і описується приростом кута між нормалями у(х,) і п(х,) до поверхні Землі ду і 


еліпсоїда (2), що перетинаються в точках х, і х або, по-суті, значеннями напрямних 


ічі? 
косинусів соб(п.х; ). Величини напрямних косинусів в точках рельєфу за профілем не визначають 
через надмірну складність таких вимірювань, тому гравіметричні виміри є лише приростами 
абсолютних значень сили тяжіння (аналітичною апроксимацією яких є значення МГПСТ) між 
поверхнями У) - Сх,. Відтак, під словами , значення сили тяжіння" слід розуміти ,значення МГ - 
ПСТ". 

Для тлумачення цих значень виконують редукцію сили тяжіння, виміряної на різних рівневих 
поверхнях, до єдиного рівня, зазвичай, рівня геоїда (принагідно зазначмо, що уточнення локально- 
го геоїда в межах України, що є окремою проблемою, теж призведе до уточнення виконуваних ре- 
дукцій). Під редукцією значень МГПСТ будемо розуміти таку задачу: визначити за заданими на 
поверхні ду значеннями «(х), хєдбу її значення «(2) в деякій точці 2 поза поверхнею ду. 

Зауваження 1. Якщо функція «(х) задана лише на частині поверхні бу (звична справа в гео- 
фізичній практиці), незалежно від того, в яку область вона продовжується, задача редукування 
хх) є некоректною |Черньій, 1992. 

В ситуації, коли для аналізу графіка за його графічним зображенням не вистачає роздільної 
здатності, застосовують його диференціальні властивості. По аналогії, для наближеного опису зна- 
чень хх) в глобальній області на додаток до їх ,,графіка" у вигляді ізоповерхонь рівня задіємо ди- 
ференціальні властивості цих поверхонь, що характеризуються відношенням (1). Для відновлення 
з гарантованою точністю поля сили тяжіння за Його гармонічним наближенням «(х) слід вирішу- 
вати відповідну граничну задачу в області малої міри. Сучасні методи |Интерпретация..., 1992; 
Гравиразведка..., 1990; Балк, 1989; Старостенко, 1978; Булах и др., 1976| вирішують такі задачі з 
точністю до деякої сталої інтегрування, що залежить від форми 1 розмірів області розв'язку. При 
відновленні значень сили тяжіння в глобальній області за її значеннями в локальних підобластях 
отримані розв'язки неможливо ,зшити", як ковдру з клаптів, через відсутність прийнятних крите- 


ріїв склеювання локальних гармонічних наближень". 

В |Алексидзе, 19585| її зведено до розв'язання внутрішньої задачі Діріхле для рівняння Лапла- 
са в локальній області; в |Черньгй, 1992) доведено коректність її постановки на множинах витоко- 
подібно зображених функцій у вигляді потенціалів простого шару, поширених на опуклій поверх- 
ні Ляпунова, і, в данину заслугам попередника у дослідженні теми, названо задачу його іменем. 

Задача ДАлексідзе у постановці Чорного |Дубовенко, 2009а| означає: знайти функцію 


Ух), х є ду, яка задовольняє в необмеженій замкнутій області у - у бу рівнянню Лапласа 


АУУ(х) «0, х є у", якщо в будь-якій точці межі ду області і на нескінченно віддаленій точці вона 


3 2 
задовольняє умовам: У атоми. 470), хе ду, М(х)-» 0 при |х|-» со, де я(х) - задана неперерв- 
Каї Хі 


на функція, у - обмежена область мас Землі, у - її необмежене доповнення без мас, ду -- ляпу- 
новська межа у |і у", така що у - КСЛу ,у зу (ду (фізична поверхня Землі). 

Вважаємо, що густина с(Є) є інтегровною функцією |Картвелишвили, 1983|, хоча вона може 
бути практично більш гладкою. Питання розв'язності 1 властивостей задачі висвітлені в |Дубовен- 
ко, 2009а; Дубовенко, 20096; Дубовенко, 2010), а ми визначимо кількісний вимір поняття , область 
малої міри", в якій справедлива властивість гармонічності функції «(х) : 

49. Оцінка міри області. Детальна оцінка міри області означає вияснення зв'язку геометрії 
області не лише з точністю вимірів і величиною аномалій, а і з геометрією еквіпотенціальних по- 
верхонь 1 точністю гармонічної апроксимації аномалій. Спробуємо дати її простіше, ніж у |Чер- 
ньй, Якимчик, 1999). Щоб оцінити по черзі вклад кожної з цих величин, згадаймо опис значень 
сили тяжіння в глобальній області |Черньгй, 1982). Класична аномалія загалом визначається як 


різниця значень спостереженого і нормального полів |Веселов, Сагитов, 196581: 


яз |в о іконі? 2 СО, з) 


в якій у(х) - напруженість нормального поля сили тяжіння, що визначається виразом 


жа У стаю Р - |агад Ц(х) со 0 - (тай Ус) -|згад т(х))сов 9 - с(х)сов 9 - що ; 


І 








де т(х)- У(х)- Й (х) - збурювальний потенціал, як різниці реального і нормального потенціалів. 
Останній визначено у близькій області у, простої геометрії (близькість областей у, і у означає 
близькість меж су, ібу , описану мірою відхилення поверхні стиснутого сфероїда су, від земно- 


го рельєфу Су), з такими ж масами, як в у , розподіленими ,, правильно " у певнім сенсі |Картве- 


СТ(х)  СУУ(х) сп(х) М ОСС 





лишвили, 1983). Істинна аномалія сили тяжіння г(х) - ш нормаль- 
бу(х) | Сп(х) ду(х)  бу(ю 
на похідна збурювального потенціалу Т(х) - виходячи з (3), дорівнює 
. 2 З(Х 
годо вд - уд -28(д5іп! 9, 4) 


де З(х) - куг між нормалями п(Х) і хх) у точці х (відхилення виска), СО05 З(Х) -(п(0, МОЮ) і 
09 - нормальна сила тяжіння від відомого розподілу мас усередині сфероїда гу, . Це визначення 


аномалій відрізняється від класичного (3), в якому аномалії є відхиленнями реального поля сили 
тяжіння Землі від його осереднення - апроксимації поля сфероїда з невідомим розподілом мас 
усередині нього. Відсутність інформації про розподіл нормальних мас ускладнює інтерпретацію 


З Точність продовжених як гармонічні функції значень сили тяжіння істотно залежить від розмірів області, в якій ставиться гранич- 
на задача, тому для отримання результату з заданою точністю нерідко доводиться відповідно зменшувати розміри області. При ре- 
дукції в глобальну область її розбивають на локальні підобласті, вирішуючи відповідну задачу для кожної з них, виходячи з місце- 
вого , нульового рівня" продовження і місцевої точності розв'язку. Ця точність залежить ще й від точності граничних даних, розмі- 
рів і форми локальної підобласті, а різниця рівнів, від яких обчислюють аномалії, може сягати кількасот мГал залежно від величини 
ундуляцій рельєфу і напруженості поля сили тяжіння. Отже, продовження в глобальну область не дозволяє обчислити в ній анома- 
лії від єдиного рівня. 


гравіаномалій. 

До речі, значення величини нормального поля можна отримати і по-іншому (немає в |Чернь(ій, 
19821), визначаючи косинус кута між нормалями до поверхні Землі і до поверхні рівня в поточній 
точці у вигляді 


со5д- (з) «У сов ідо мл) иа) ', 1 ато) і а ето 








РАН 2) 2() ду | е(х ду 





Рис. 2. Геометричне тлумачення аномалій. 
от(х) у з 

звідки у(х) а з (сов 30. Прикметно, що зображення аномалії (3) не є однозначним; його 
у 


0 дт(х) 
п 


альтернативна подача да (х) з М) 8іп гл аналогічна за точністю, але має інший напрям 


нормалей у і п. Оскільки вектор «(х) має унікальний нефіксований напрям в кожній точці 
поверхні бу, то обчислення за (3) майже завжди відрізняються від істинного поля. Геометрична 
інтерпретація істинної аномалії очевидна з рис. 2, де СЕ -- проекція Да (х) на нормаль у, СР - 
відхилення класичної аномалії ба (х) від проекції СЕ. Отже, є-ОВ-ОРр, у-окЕ, 


Де з ВЕ-ОБ-ОВ- | у-8|, дез рЕ- Й - | 5, . Істинна аномалія має ясне аналітичне зображення 


де) -|в)- Сг) «а |иО" - до) У сов. 


Виражаючи величину хх) через вираз (3) 1 підставляючи її в останню рівність, отримаємо: 














ДЕ НК Б З (ак 
лях) - ДУ(х), |зіп' -- я зі - з (3) 
и ЕС 
Навіть якщо класична аномалія дорівнює нулю, істинна відрізняється від неї на величину 
М 2 
Де (2) - 3) віп го (6) 


За відносної зміни нормальної сили тяжіння від полюса до екватора в 3000 мГал і варіації від- 
хилення виска 0"--40" така різниця (6) виливається у величину до 200 мГал (табл. 1), в якій цілком 
може ,уміститись" аномалія від глибинних неоднорідностей. Це підтверджено тим, що експери- 
ментально отримані для деяких районів Землі відхилення виска |Іегсі еї а!., 1979| не можна пояс- 
нити лише видимою нерівністю рельєфу, доводиться постулювати їх зв'язок з аномальною будо- 
вою кори і мантії. 


Табл. 1. Залежність розмірів | локальної області (км) від похибки гармонічної апроксимації Є, і амплітуд Т, ано- 


малій сили тяжіння (мГал). 
































Похибка Амплітуда аномалій Т,, мГал 

ві 10 |20 30 40 50 100 | 150 | 300 | 500 
0.1 42 122 14 

0.2 84 (142 28 22 16 

0.5 212 1106 |70 54 42 22 

1 424 1212 1142 1106 | 86 42 28 10 

2 142 | 170 84 56 28 10 
5 212 1142 |70 42 
10 424 |284 | 142 84 
20 950 | 568 |284 | 170 






































Прикметно, що розміри локальних областей чітко визначені і не надто великі; за межами окре- 
слених ними зон закінчується сфера впливу класичних трансформацій і починається сфера пливу 
нелінійних наближень, отриманих завдяки розв'язанню задачі Алексідзе. Звичайно, наведені в 
табл. 1 результати приблизні, оскільки врахована лише так звана операторна похибка є,, а вплив 


на оцінки здійснює і похибка вимірювань є,, і похибка округлень є,, які мають кумулятивний 
ефект. 

Відзначено |Алексидзе, 1985|, що якби аномалія ба (х) була високоточною, то перерахунок у 
зовнішньому просторі був би неможливим через істотну залежність рівняння сили тяжіння від 
швидкості зміни кутів між напрямом аномалії і координатними осями, тобто величини а(х) ГЧер- 
ньй, 1970). Фіксувати швидкість такого обертання неможливо, а нехтування цим відхиленням іс- 
тотно спотворює перерахунок (на величину (б6)). В такій ситуації доречне припущення Юнга, оче- 
видне з рис. 2: класична аномалія з достатньою точністю, що не перевищує величини 
СР «ОК-ОС або є, -2 х(х) віп? 0/2 чисельно співпадає з проекцією вектора аномалії на на- 


Т(х) 


прям вектора Ух) , тобто до (х) з цю - Є, (табл. 1). 
у 


Розглядаючи аномалії дах) з точністю до є,, ми відкидаємо величини порядку 1- со5(п, у), 


3 
де  1- со5(п, у) -2віп! - У сов(л, Хі; Їсоз(х, , п) Е5 соз(х, ; у)| б У соз(у, Хі сов, 5 у) М соз(х, 5 п)|. 
ке СІ 


Якщо нормальний потенціал обрано вдало, то в замкнутій необмеженій області ду - в? ху 


Т(х) 


99 


завжди виконується умова б б тах, (о) ««тах 
хедбу хебу 


, Тому можна знехтувати вказаною величи- 








ною. 

На основі наведених даних приходим до висновку, що з точністю до порядку 0.001 « є, 2 0.05 
мГал аномалії сили тяжіння до (х), х є бу є значеннями модуля проекції градієнта збурювального 
потенціалу т(х) на напрям нормалі у до еквіпотенціальної поверхні нормального поля сили тя- 
жіння, тобто 








х 
чи ба (х) - |огад т(х) «сову, х, | - бе (х) кара (7) 
хебу у 
У прийнятій нами системі координат напрям нормалі співпадає з віссю 7, тому справедливим 
ОТ(Р бо й 
буде відношення да (х)- гир) - єр у будь-якій точці Р(х) деякого околу у,. Якщо аномалії ви- 
с 


2т(2) 


дх. 


1 


значити з точністю є,, а складники градієнта по горизонталі «є, » 0, не перевищують 
ізі,2 


заданої константи, то за умови, що вклади кожного з напрямних косинусів со5(х,, у), і-12 1 








1-- со8(хз, у) однакові, легко знайти мажоранту для виразу (7): 


Є Ж 8 


сов(ху)я 1, гої; (8) 





со5(х, У) -Д| є дб, |сов(х, м є 





3є, 


У визначеній таким чином області аномалії сили тяжіння да (х) з точністю є, - є; можна вва- 


жати гармонічними функціями. З геометричної точки зору обмеження (8) є умовою обертання нор- 
малі (Р) навколо осі 2, спрямованій по нормалі у(х) до еквіпотенціальної поверхні (/ (х)- С нор- 


мального поля в декартовій системі координат з центром в точці х. При русі навколо нормалі у(х) 


локальна область І-Всо50 





екваторіальна площина 


Рис. 3. Геометричне тлумачення форми і розміру локальної області гармонічності гравітаційних аномалій. 
нормаль у(Р) описує в першому наближенні поверхню кругового конуса. 

Якщо вершину цього конуса помістити в точку перетину нормалі у(х) з еквіпотенціальною 
площиною еліпсоїда (рис. 3), то локальна область буде обмежена бічною поверхнею конуса, який 


місить нескінченно віддалену точку по осі 2, та земним рельєфом. З рис. 3 очевидно, що радіус 


ве 


вирізаної на поверхні Землі області бу, не перевищує величини І - Кд, де да 1, К - радіус 


З, 
Землі в поточній точці Р. Згідно правил диференціальної геометрії, конус, зорієнтований як на 
рис. 3, матиме вершину в центрі середньої кривизни нормальних перетинів еквіпотенціальної 
поверхні (/ (х)- С в точці 0, тому у виразі / - Кб радіус К характеризуватиме середній радіус 
кривизни нормальних перетинів поверхні (/ С 
Виходячи з (7) 1 (8), приходимо до висновку, що діаметр основи локальної області не переви- 


щує величини 4(х) «21 - 2К(є, ч є, ) зтах що 
хесу Х з 





, тобто діаметр плоского околу ду, області вимі- 
рювань при відомих точності вимірювань є, та апроксимації є, (ці апріорі відомі величини визна- 


. З з ТІ х 
чаються методикою спостережнь) тим більший, чим менша амплітуда пах дт(2) 


спостереженої 
хеду дх. 


аномалії, а при фіксованій амплітуді росте разом зі зібльшенням похибки вимірювань. Нескладні 
розрахунки показують, що для того, щоб описати аномалію сили тяжіння амплітудою 30 мГал з 
точністю 0.1-0.5 мГал, слід не виходити за межі зони а -8-35 км. 

Для точної оцінки розмірів локальної області необхідно точно знати аномальні висоти геоїда. 
Для уточнення фігури Землі за формулою Брунса необхідно з високою точністю визначити збурю- 
вальний потенціал Т(Х); це можна зробити, наприклад, за алгоритмом (|Черньгй, Якимчик, 1999) у 
вигляді розв'язання послідовності задач Неймана із спеціальними граничними умовами. У першо- 
му наближенні можна скористатись прийомом, наведеним в | ДДубовенко, 2008а, 6. 

Зауваження 2. Якщо для перетворення аномалій да (х) використовують конструкції з нескін- 


ченними межами по осях х,, і- 1,2 (типу інтегралів Пуассона), то розміри локального околу уз- 


довж осі г різко скорочуються до 0.21. 
59, Похибка трансформацій. Якщо нормальний сфероїд 1 розподіл мас у ньому обрані вдало, 
відхилення виска у районах зі ,спокійним" рельєфом не перевищує 60", і оцінка зверху останньо- 


го доданка у (4) не більша 2о(х)віп" пи «0,05 мГал, хєсу. Зв'язок геометрії локальної області з 


кривизною еквіпотенціальних поверхонь знайдемо через різницю у околі точки хє у" функцій 

и(у) 1 у(у) продовження аномалій (у), ує сгу(х), розв'язавши відповідні задачі Діріхле для 
рівнянь аномалій сили тяжіння 1 Лапласа |Черньгй, Якимчик, 1999): 

Ди(у) -а"(у)и(у)є 0, уєу(х)у",и()-тбу), уєду(ю, (9) 

ДУ(у)є0, уєУ(Д0У", У) є (У), УЄ 2У(2. (10) 

Задачі (9) і (10) співпадають при а"(у) -0, у є У(.) (вектори одиничних нормалей п(х) у дові- 

льних точках Х паралельні, наближення до такої ідеалізації є в областях малої міри - в них напря- 

мок , повного" вектора сили тяжіння постійний, а відносні флуктуації сили тяжіння від середнього 


її значення невеликі). В іншому разі оцінка точності трансформацій гравіаномалій як гармонічних 
функцій має вигляд 


з а"(при(п) 
- РУ брус | | 
є з тах|є 15, - шах, Т У-л і 7 





12) 


з: 
с 


т а" 
г з тах |(УТО(владу бо) «7-Й | 


(а, Б - велика і мала півосі сфероїда С, ). Похибки продовжених аномалій сили тяжіння внаслі- 


док заміни субгармонічної функції гармонічною функцією залежить від кривизни еквіпотенціаль- 
них поверхонь поля, величини аномалій і міри локальної області, у яку здійснюють продовження. 
Похибка трансформації аномалії г(у), ує сУ(х) за (11) для обчислення параметрів будови ша- 


ру Землі товщиною 200 км в обмеженій області У(х) (усічений консферичними поверхнями круг- 
лий конус із вершиною у центрі Землі, основи якого віддалені від центра на К і К,, верхня основа 





співпадає з поверхнею ду, і має діаметр 0--2200 км, х 4200 - К - 6371,1 км - середній радіус Зе- 


млі) дорівнює 





2 «2143 сові зро -дгсовфаг Я 
То 


о ГОЛ -совф нг) ро -дгсовфаг -1 ро -дрсовфаг -со8Ф 


со8 фіп і 
1 -со5Ф ро -дгсовфагі І 1--42Й - со8 фФ) - со5Ф 


- амплітуда аномалії, г-К/К. 














де г, - |тах г(у)- посі у) 


уєу(О) уєУ(Д 








Для локальної області У(х) є/ т, з 0,5:107 --0,3-10" і за відомої амплітуди аномалії легко знай- 
ти похибку відхилень субгармонічної функції від гармонічної, але амплітуди аномалій невідомі. 
Нехай у локальних областях типу у(х) т7,210000 мГал, тоді 5 4 є 2290 мГал. 

Похибка продовжених аномалій г(х), хє су, як гармонічних функцій в область у(х) пу", За 
умови, що амплітуди продовжених аномалій 7, 2150 мГал, а локальні області У(х) - усічені кону- 
си (0930! « ф «109; 0 « В « 3602; 637114 р 6571, км), нижні основи яких співпадають із земним 
рельєфом, - за (11) складає 0,07 « є 46,3 мГал. За (12), якщо діаметри нижніх основ сімейства усі- 
чених конусів 110«|бх|«2200 км сфероїда Фу: а-6378,2 км, Б -6356)9 км, то 3,7 4 є 4 73,4 
мГал. Оцінки (11) і (12) різні в області УС)п у" через те, що в них замість кривизни аЧт) граві- 
поля в точках 7 є У(х) (залежної від особливостей регіонального поля і поточної аномалії) вжиті її 
наближення а; (п) йо (п) у вигляді концентричних сфер і конфокальних еліпсоїдів. 

Область малої міри у(х) перебуває поблизу локальних і на достатній віддалі від регіональних 


особливостей поля, тому птаха; (17) « таха"(п), і- 0,1, тобто, обидві оцінки похибки занижені. 
леу() леу() 


Похибка (11) аг (п) З а (п) «0,5-10"" істотно менша оцінки (12) і аз(п) « У а: (а,0,0) «0,22:107: 





у першій формулі вжито наближення а? (п)), а у другій - ча; (п). 

Оцінки похибок трансформації в кожній локальній області У(9 і У(у), навіть перетнутих між 
собою, не дають змоги оцінити загальну похибку трансформації в об'єднанні у(х) 7 У(У), 
ху у є К?). Якщо у точці х відхилення виска 4(х)-60", а в точці у-Ху)-0, то ,істинна" 
аномалія дбе(х) в точці х, у вигляді модуля різниці градієнтів реального й нормального 
потенціалів, дорівнює 





бе(х) з (0) - у02| -- НУ віп? м -4у(0г(9 віп! | зт), 


ауточціу бе(у) - т(у). В переході з локальної області У(х) до області у(у) істинна аномалія від- 


З(х) 


різнятиметься від класичної на 62(х) - 2/(х) зар ж 285 мГал, навіть якщо класична аномалія 


то)с0 в точках х і у. Трансформація поля г(/)) (Ди(7)) - 0) в кожній з областей У(х) і уУ(у) 
при г(х)-т(у)-0 і х у як гармонічних функцій у(п/) не дозволяє виявити істинну аномалію, 


яка сягає величин, порівняних зі значеннями трансформованих аномалій навіть у районах зі спо- 
кійним рельєфом. 

Отже, перетворення поля аномалій сили тяжіння в різних локальних областях за правилами 
трансформації гармонічних функцій не дає можливості виявити істинну аномалію, яка навіть в ра- 
йонах із незбуреним рельєфом може бути р. ітик ною з редукованою. Саме тому область за- 
стосування методів інтерпретації аномалій сили тяжіння, розвинутих на базі теорії гармонічних 
функцій, обмежується рамками розвідувальної геофізики; їх не можна без додаткових обмежень і 
удосконалень застосувати для вивчення глибинної регіональної структури Землі. У зв'язку з цим 
риторичним є питання, враховувати чи ні природу аномалій при глибинному вивченні Землі і не 
має альтернативи задача Алексідзе | Дубовенко, 20094|. У світлі сказаного вона є однією з найакту- 
альніших проблем гравіметрії в ХХІ ст., 1 повністю вкладається в схему розвитку геофізики (аналі- 
тичні апроксимації середовища і поля, нелінійні задачі, дискретні постановки, великорозмірні сис- 
теми рівнянь) Страхов и др.., 2001. 

Автор висловлює глибоку вдячність академіку НАН України В.І. Старостенку за цінні заува- 
ження. 
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